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We consider the antiferromagnetic spin- 1/2 Heisenberg model on a two-dimensional bipartite 
quasiperiodic tiling. The broken symmetry ground state in this model is inhomogeneous, but never- 
theless bears interesting similarities with that of the square lattice antiferromagnet. An approximate 
block spin renormalization scheme developed first for the square lattice is generalized here to the 
quasiperiodic case. The ground state energy and local staggered magnetizations for this system are 
calculated, and compared with the results of a Quantum Monte Carlo calculation for the tiling. 
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I. THE EXPERIMENTAL AND THEORETICAL 
BACKGROUND 

Magnetism in quasicrystals can be very complex, due 
to the extreme sensitivity to structural details, in such 
systems, of local moment amplitudes as well as of the 
interactions. A considerable simplification of the problem 
is however possible for the recently studied rare-earth 
based quasiperiodic alloy ZnMgHo Q. The rare-earth 
based magnetic alloys represents a conceptually simpler 
system than the transition metal alloy quasicrystals that 
were initially the object of experimental studies, since the 
magnetic moments are associated with f-orbitals, and can 
be assumed in the first approximation to be independent 
of the local itinerant-electron density of states. 

This is to be contrasted with the earliest magnetic qua- 
sicrystals of the AIMn family, where the itinerant mag- 
netic moments on the Mn atoms depend sensitively on de- 
tailed structural features due to the d-orbital hybridiza- 
tion (see the review by Hippert et al in 0). To add 
to the difficulties the early alloys were metastable qua- 
sicrystals of inferior structural quality so that the role of 
disorder had to be considered in addition to the intrinsic 
behavior. Experimental results indicated a wide distribu- 
tion of effective moments on the Mn atoms, as well as of 
the interactions between these, leading to a large number 
of unknown parameters in the phenomenological models 
describing such systems. From a theoretical viewpoint, 
therefore, the rare earth system is clearly far simpler. 

ZnMgHo was shown to undergo a magnetic transition 
into a magnetic state characterized by short range an- 
tiferromagnetic correlations with quasiperiodic modula- 
tion The experimental results lead naturally to the 
question of what properties one expects for the ground 
state of a quasicrystal with short range antiferromagnetic 
interactions. An acceptable starting point for models 
of such systems could be, as for crystalline compounds, 
a Hamiltonian with short range antiferromagnetic cou- 
plings between pairs of identical spins, H = ^ Jy Si.Sj. 

Fig.l shows the results of a recent Monte Carlo study 




FIG. 1: Inhomogeneous ground state structure on the tiling. 
The circles have sizes that depend on the strength of the local 
order parameter 



of a two-dimensional model of quantum spins on a 
quasiperiodic tiling The circles on the vertices have 
radii that depend on the value of the local staggered mo- 
ment, a quantity that we will define further below. The 
tiling considered is the eight-fold symmetric octagonal 
(Ammann-Beenker) tiling, in which sites can have six 
possible values of coordination number z. Sites were 
occupied by S = | spins, with uniform interactions 
Jij = J > along the edges of the tiling. The system is 
bipartite, meaning that every spin belongs to one of two 
sublattices and interactions couple only spins of differ- 
ent sublattices. Analogously to the spin ^ square lattice 
antiferromagnet, which is now believed to have a ground 
state with long range order, we expect that this quasiperi- 
odic system, too, has a broken symmetry ground state 
with long range order. Classically, the ground state cor- 
responds to having oppositely directed sublattice magne- 
tizations, with no frustration, in the sense that all bonds 
can be "satisfied" simultaneously. In the quantum case, 
the ground state will correspond to zero total spin since 
for the octagonal tiling, the two sublattices are equiva- 
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lent. 

The inhomogeneous structure of the ground state seen 
in Fig.l is a reflection of the environment-dependence of 
the quantum fluctuations around the Neel state in the 
quasicrystal. There is for the moment no spin wave ex- 
pansion that would allow, as in a periodic solid, to calcu- 
late the distribution of staggered moments and explain 
the QMC results. In fact, as we will see, a real space ap- 
proach seems more appropriate for quasiperiodic tilings, 
and many such calculations exist for the case of one di- 
mension. 

One-dimensional models to study the behavior of 
quantum spins on quasiperiodic chains been considered 
by several authors. Quantum spin chains have been an- 
alyzed using renormalization schemes 0, based on the 
inflation symmetry of these chains. Using a mapping to 
fermionic models and techniques of bosonization, |6j it 
is possible to obtain interesting results concerning global 
properties such as the magnetization as a function of ex- 
ternal field, and the spectral gaps for a variety of dif- 
ferent quasiperiodic sequences. However, real space in- 
formation such as the distribution of the local quantities 
mf oc i ~ has not so far been calculated. 

For two dimensional structures, real space configura- 
tions have been studied for models with classical spins. 
Here, the ground state is nontrivial only when the model 
includes frustration. In Q Godreche et al introduced 
a renormalization scheme on the Penrose tiling for a 
Heisenberg exchange model with competing antiferro- 
magnetic interactions, and were thus able to obtain a 
phase diagram consisting of a variety of ordered phases. 
The real space spin configurations were calculated numer- 
ically in [8j for classical spins interacting via long-ranged 
dipolar interactions, and a complex magnetization dis- 
tribution with overlapping decagonal rings reflecting the 
underlying Penrose tiling was found. A quasiperiodic 
magnetic state with a heirarchical ringed structure was 
found, as well, in a different context: that of itinerant 
magnetism due to interacting electrons 0. 

With this background, we return to the problem of 
quantum Heisenberg spins with nearest neighbor antifcr- 
romagnetic interactions. Ref. Q presented local stag- 
gered order parameters similar to the quantities ro; OCj j 
above, calculated for individual sites using the expecta- 
tion values for local spin-spin correlations. One sees in 
Fig. 1 that sites of the same z have similar local order 
parameter amplitudes. An explanation of this behav- 
ior was given by considering isolated star shaped clusters 
called Heisenberg stars in [2j . This provided a qualitative 
understanding of the decrease of local staggered magneti- 
zations as a function of z, but for a more quantitative fit 
to the QMC results, it is necessary to go beyond the iso- 
lated cluster approximation, and take into account longer 
range correlations. This can be done in a renormalization 
group (RG) calculation that uses an important symme- 
try of the tiling, namely invariance under discrete scale 



transformations called inflations. This renormalization 
group is a generalization of the calculation of Sierra and 
Martin-Delgado for the square lattice 0| , where the au- 
thors considered star-shaped block spins formed by a cen- 
tral spin and its four nearest neighbors. In their calcu- 
lation, block spins formed from these five-spin clusters 
are shown to interact via an effective Heisenberg anti- 
ferromagnetic interaction on a bigger square 
lattice. The effective spin values scale to infinity, i.e. the 
classical limit, under renormalization. Their model for 
a translationally invariant system can, as we will see, be 
adapted to our quasiperiodic case under certain approx- 
imations. We thus calculate not only the global ground 
state energy as was done for the square lattice, but also 
the distribution of local order parameters. We will dis- 
cuss the method, which has been briefly reported in , 
in some detail in the present paper. 

We begin with an introduction to the quasiperiodic 
tiling and the spin Hamiltonian in the next two sections. 
The RG scheme is described in the fourth section. Re- 
sults and discussions are presented in sections five and 
six. 



II. REVIEW OF GEOMETRICAL ASPECTS 




FIG. 2: A portion of the octagonal tiling showing the six 
different nearest neighbor environments A,B,...,F 



1. Some general remarks 

The octagonal tiling ^1 shown in Fig^Jcan be thought 
of as the equivalent of the square lattice for quasiperi- 
odic systems. It has therefore been frequently used for 
analytical and numerical investigations of the effects of 
quasiperiodic modulations in two dimensions. Spectral 
properties of electro ns Il3| , transport properties Il4j , vi- 
brational properties [15j and magnetic properties [2j have 
thus been studied for discrete models defined on the oc- 
tagonal tiling. The tiling is built from two kinds of tiles, 
squares and 45° rhombuses. These two types of tiles can 



fill the two-dimensional plane in an aperiodic way, as 
Penrose first showed for the five-fold tiling named after 
him dHH. 

Although there is no translational invariance in a 
quasiperiodic tiling, any given tile arrangement of tiles 
reoccurs all over the tiling with a certain frequency of 
re-occurrence - or, alternatively viewed, there exists a 
mean distance of separation between such identical do- 
mains. This is referred to as the repetitivity property 
of quasiperiodic tilings, and is very different from the 
situation in a disordered medium (where the expected 
distance in which one expects to find a second region 
identical to the first increases exponentially with the size 
of the region) . Similarly, the property of symmetry un- 
der rotations for these tilings differs from that in crystals, 
for which the new and the old structures coincide exactly. 
For the quasicrystal, the equivalence of the new and old 
tilings holds in the "weak" sense, namely, any finite re- 
gion of the new tiling after rotation will be identical to 
finite regions of the old one. 

Such aperiodic structures can be built using "match- 
ing rules". These are local rules that determine if and 
how two tiles can be laid side by side (see Ch.l of (lif). 
Alternatively, tilings such as the Penrose and octagonal 
tilings could be generated by a projection method down 
from a higher dimensional periodic structure 19]. Such 
an approach can give either a deterministic, perfectly or- 
dered tiling, or a random one where tiles are assembled 
subject only to the constraint that they should fill space 
without overlapping [2(Tj . Random tilings are of great 
theoretical interest, but we are here interested in deter- 
ministic tilings, which have the important property of 
invariance under inflation/deflations, or discrete scale in- 
variance. This symmetry is illustrated in Fig. |3|and will 
be described in more detail in the next section. It is this 
property that is responsible for the characteristic singu- 
lar electronic and magnetic properties of such tilings and 
it was first pointed out in the Penrose tiling, which is in- 
variant under a replacement of tiles by r-fold bigger tiles, 
where t = (\/5 + l)/2 [l7|. One can define geometrical 
inflation rules for, among others, the Fibonacci chain in 
one dimension, the octagonal tiling in two dimensions, 
and the icosahedral tiling in three dimensions. 

The renormalization approach is a natural one for such 
geometrically self-similar quasiperiodic tilings , and this 
structural property has been exploited in order to es- 
tablish recurrence relations for parameters occurring in 
discrete spin models, electron hopping models, etc, as 
mentioned before for the one-dimensional case, but also 
for some two-dimensional models 0, EI , where analyti- 
cal methods remain hard to implement. As noted in the 
introduction, our approach is inspired by the renormal- 
ization calculation of Sierra and Martin-Delgado 0] for 
the square lattice. 

Some principal properties of the octagonal tiling that 
are used in the RG calculation are reviewed in the 
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FIG. 3: Portion of original (black) tiling, showing sites of the 
a class (black dots) which become sites of the new inflated 
(grey) tiling 



next section, without demonstration. (For those inter- 
ested, Appendix A contains some additional details on 
how to obtain a quasiperiodic structure, and how infla- 
tions/deflations are described in the framework of the 
projection method. Although not strictly necessary to 
understand the calculations presented below, an under- 
standing of the geometrical properties of the tiling is im- 
portant for those wishing to improve this approximate 
RG scheme and extend it to other models. For more 
details, the reader is referred to reviews in [22Ll23| ). 



2. The six local environments 




FIG. 4: The six different nearest neighbor environments of 
the octagonal tiling 

The six nearest neighbor configurations, correspond- 
ing to coordination numbers z — 8, 7,..., 3 are labeled 
A,B,...,F as shown in Fig. 2. Fig. Qlshows these environ- 
ments separately. In an infinite tiling, each of these types 
of site occurs with a well-defined frequency fi , where (see 
Appendix A) 

I A = A~ 4 ; fs = A~ 5 ; f c = 2A~ 4 ; 
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foi = A" 3 = f D2 J E = 2A- 2 ; fp = A" 2 (1) 

with A = 1 + \/2. One distinguishes between two kinds 
of D sites as explained in the next section. It can be 
checked using the above frequencies that the average site 
coordination number on the octagonal tiling is exactly 
four. 



3. The inflation transformation 

Inflation proceeds as follows for the octagonal tiling: 
one starts with a tiling composed of tiles of a given ini- 
tial edge length (we will assume this is equal to 1) and 
one reconnects a precisely determined subset of vertices 
so as to obtain a new tiling of the same type as the old, 
i.e. having the same set of local geometries, except for 
an overall scale change by a numerical factor A = 1 + \/2 
( Fig. 4). The sites shown as black dots on the origi- 
nal tiling belong in the a class: A,B,C and half the D 
(called Di) sites. These become the sites of the new big- 
ger tiling, while the remaining (/3) sites drop out. Note 
that there are two varieties of five-fold sites, D\ and D 2 , 
which belong to the a and (3 classes respectively. On the 
octagonal tiling, they always occur in pairs. Appendix A 
shows how the two classes of D sites can be distinguished 
in terms of their perpendicular space coordinates |2^ . 

Under inflation, the density of sites is reduced to 1/ A 2 
of its initial value. The sites that remain acquire new 
values of the site coordination numbers z' < z. The table 
below lists the initial and final values of coordination 
number for each of the a class sites (note that there are 
four different subcategories for the A sites - see Appendix 
A for more on the properties of these subcategories). 



initial site 


final site 


(z value ) 


(z value) 


A (8) 


-> A,B,C or Di (8,7,6,5) 


B(7) 


D 2 (5) 


G(6) 


E (4) 


Di (5) 


F(3) 



Table 1. List of a sites and their transformations under infla- 
tion 



4. Nearest neighbors of a sites 

For the four types of a sites, the table below lists the 
nearest neighbors (nn) in terms of the type of site and 
the number of sites of that type. This information will 
be useful in determining the final block spin value at the 
central site, as we will explain in sec. III. 



a site 


nn spin type(number) 


A 


F (eight) 


B 


F (five), E (two) 


C 


F (two), E (four) 




D 2 (one),E (four) 



Table 2. The a sites and their nearest neighbor 
environments 

Table 2, in conjunction with Table I allows one to de- 
duce how blocks are organized in the tiling. An A site 
which transforms to an A site after inflation corresponds, 
on the original tiling to an A block surrounded by eight 
D\ blocks. 

III. THE SPIN HAMILTONIAN 

We consider onsite spins Si ( i=l,N) where all spins 
have spin ~, with the Hamiltonian H(N, {Si}, { }) 

H = JySf.Sj (2) 

(id) 

where denotes a pair of spins Si and Sj linked by 
an edge, and = J > for such a pair. This system 
is bipartite with two identical (to be understood in the 
weak sense) subtilings (as on the square lattice) . 

1. Finite spin clusters 

A sites are surrounded by eight F sites. If one isolated 
one such cluster of 8+1 spins, the lowest energy state for 
classical spins is the one with the eight peripheral spins 
antiparallel to the central spin. In the quantum case, the 
ground state of the cluster is rotationally invariant, and 
corresponds to the total cluster spin value S to t = 7/2. 
The other a sites correspond to total cluster spin values 
in the ground state of S to t = {zb — l)/2 = 3 around a B 
site, and so on. The four clusters are shown in the left 
hand side series of FigJSJ 

Clusters of each type can be defined on larger and 
larger length scales, by using the inflation rules already 
outlined to determine the new A,B,C and D\ sites after 
inflation. FigJSJ shows the four a clusters on the next 
largest length scale on the right hand side series. Here, 
block spin centers are shown with big black dots, while 
the sites corresponding to the (3 sites are indicated by 
smaller dots. On a yet bigger length scale, FigElshows a 
"second generation A site" , namely, a site that remains 
of A type after two inflations, along with all the sites 
belonging to the cluster before the two decimations. 

The a clusters on all length scales are the building 
blocks for the renormalization scheme that follows. 





FIG. 5: The a site clusters defined on the original (left) and 
once inflated (right) tilings 



IV. THE RENORMALIZATION 
TRANSFORMATION 



The renormalization calculation is a generalization to 
an aperiodic system of the one used for the square lattice 
by Sierra and Martin-Delgado [lfj. We review briefly 
the steps of their calculation before showing how they 
are modified in the quasiperiodic case. 



FIG. 6: Second generation A cluster 



1. RG on the square lattice 

We consider the nearest neighbor Heisenberg antifer- 
romagnet described by Eq. 2 with spin ^ on the vertices 
and the initial coupling J along the edges of the squares 
(of side a = 1). Fig[7| shows the five-spin blocks en- 
closed by circles. The four couplings inside each block 
are shown outlined by thick grey lines. As one sees, the 
block spins form a new rotated square lattice of side 
(FigQ). Each of the blocks can be diagonalized exactly. 
With every step of RG, only the lowest energy states of 
the blocks are retained to form the basis for the effective 
Hamiltonian. To and Tq denote the operators describ- 
ing the transformations from the original Hamiltonian 
(acting in the complete Hilbert space) to the effective 
Hamiltonian (acting in the reduced Hilbert space). For a 
single block, the lowest energy sector corresponds to spin 
|, and the ground state energy is eo = — J5(45 + 1). 
The couplings not already taken into account give rise to 
inter-block interactions, calculated by first order pertur- 
bation theory. It is easy to check that the new block spins 
will be coupled antiferromagnetically to its nearest neigh- 
bors, like the original spins. The effective Hamiltonian 
H(N, 5, J) can thus be written approximately as a sum of 
single-block contributions (the diagonal terms) and a set 
of terms involving nearest neighbor blocks (off-diagonal 
terms), and the formal expression for the transformed 
problem reads 



T^H(N, S, J)T Q = jV'e (J, S) + H'(N', S' , J') 



(3) 
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where the new Hamiltonian H' has the same form (bi- 
linear in S') as H , and N' — N/5. The effective spin of 
a block spin is S' = 35 = | . The spin renormalization 
factor relating one of the four boundary spins to the new 
block spin has been shown to be close to the classical 
value £o = Si/S' » | ( see for the exact value) . The 
interaction between two contiguous blocks is J' — 3£q J. 

Repeating the steps of renormalization, one has ulti- 
mately for the ground state energy per site an infinite 
sum as follows 

oo 

eoo = — ^5-"J w 5 w (4x3 , '5+l) (4) 

5 n=0 

where S^+V = 35^ and = 3£ 2 (S"^) J (n) . 

Under RG, the spins evolve to the classical limit, S — > oo 
indicating that in the quantum case as well one has 
a ground state with broken symmetry. The couplings 
scale to zero indicating the model is massless. Qualita- 
tively, thus, the RG gives the now accepted physics of 
the model, however, quantitatively the value obtained 
for eoo w —0.546 is not as good as that obtained by spin 
wave expansion and is about 15 % higher than that es- 
tablished by numerical calculations [25j. We will return 
to this point at the end of the paper. 
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FIG. 7: Five-spin units (surrounded by circles) on the square 
lattice. The new vo x \/5 unit cell is shown 



2. RG on the octagonal tiling 

On the octagonal tiling, it is clear that several kinds 
of block spins must be introduced. A natural choice is to 
designate the a sites as block spin centers. Figll2l shows 
the positions of the block spins (black dots) on a portion 
of the tiling. Upon inflation, the other sites will disap- 
pear, leaving only the block variables, and some residual 
interactions between them. If no new couplings are gen- 
erated, one will find an effective Hamiltonian similar to 



the old, except for the renormalized couplings which be- 
come site dependent. One can repeat the process, and 
determine if there is convergence to a fixed point. 

The simple scheme outlined above cannot be imple- 
mented without some modifications and approximations. 
The first problem arises because the connectivity of the 
tiling is such that some of the block spins overlap, that 
is, share two intermediate f3 sites in common. This is 
shown by the thick grey lines in Fig^J which indicate 
the boundary between overlapping blocks. Overlapping 
occurs between contiguous C and D\ blocks, as well as 
between contiguous D\ blocks. This overlapping occurs 
with a finite density. One can calculate this density by 
noting that the shared sites occur between any two sites 
that are a distance X 2 d s apart, where d s is the short 
diagonal of the rhombus. One finds, using the relative 
frequencies of occurrence of squares and rhombuses that 
the density of pairs is \/2/A 3 , that is, about 10% of the 
total number of pairs. 




FIG. 8: Tiling showing block centers (black dots). The grey 
lines connect pairs of sites that are shared between two blocks. 



To deal with this problem, we therefore considered two 
possible modifications of the original model, i) doubling 
the number of spins on each shared site, and considering 
each spin as being coupled to one block only, and ii) de- 
coupling the block spins by annulling one of the bonds to 
the left or the right so that spins are no longer coupled 
on both sides. The first modification leads to overesti- 
mating the total energy, the second to underestimating 
it, with respect to the original octagonal tiling. Spin 
doubling on selected sites leads to an uninteresting flow 
under renormalization, where cluster energies basically 
repeat a scaled Heisenberg star distribution at each step. 
The bond dilution scheme yields a more complicated be- 
havior of cluster energies under renormalization, and is 
the option taken up in detail in this paper. 

We note that the diluted model remains two- 
dimensional, and is not of a scale invariant fractal such as 
the Sierpinski gasket [i(J, where bonds are also deleted 
heirarchically but in a way that leads to an effective frac- 
tal dimension less than two. 

The second problem is the quasiperiodic connectivity 
between blocks which leads ultimately to an infinite num- 
ber of environments. This is dealt with by truncating 



the number of environments we choose to distinguish be- 
tween. The a sites always have the same type of nearest 
neighbors (given in Table 2), however the sites occur 
in several configurations. We will now truncate the table 
of connectivities by allowing only one type of Di ,E and 
F site, and a connectivity table as follows: 



(3 site 


nn spin type(number) 


D 2 

E 

F 


Di (one),E (two),F(two) 
a (two), F (two) 
a (one), E (two) 



Table 3. The f3 sites and the truncated set of nearest 
neighbor environments 

A. Bond dilution and the new block spins 

In this subsection we discuss the blocks that are ob- 
tained after dilution and the values of the effective block 
spin. Fig[!Jtop) shows in detail a central D\ site which 
transforms to an F site under inflation. The three neigh- 
boring block spins are shown as well, with the block 
spin sites shown by black dots. The original links are 
indicated by thin black lines, while the new effective 
links on the inflated tiling are shown by thick grey lines. 
FigOJmiddle) shows a C site transforming to an E site, 
with the same conventions used to denote block spin sites 
and new effective couplings. In this figure one sees that 
two of the block spins, corresponding to neighboring D\ 
blocks, overlap. The pair of sites shared between the two 
blocks is coupled to the left and right by a total of four 
bonds. In the bond-dilution approach, one has to set 
two of the bonds equal to zero. This can be done in one 
of two ways that treat the two blocks equitably, leaving 
each D\ block with one less bond. Finally, Fig[!|f bottom) 
shows an A site transformed under inflation to a final A 
site. In this case, the eight D\ blocks surrounding the 
center block form a ring of overlapping blocks. There are 
two ways to decouple them all by annulling eight of the 
sixteen links joining them in way that treats all the D\ 
sites equitably. Ultimately, the bond dilution results in 
an effective reduction of connectivity of C and D\ sites: 
the former have the effective z value z = 5 and the latter 
5 = 3. 

B. Spin renormalization factors 

Consider a block spin composed from a cluster of z 
spins surrounding a central spin and antiferromagnetic 
interactions. In the simplest case where all spins have 
the value S, the block has a spin of S' — (z — 1)5 in 
the ground state. The spin renormalization factors are 
taken to be equal to the classical value for simplicity, 
so that for a given block £ z = (z — The new 
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FIG. 9: Block spin centers (filled circles) showing the central 
and all peripheral blocks for three cases: (top) a z = 5, z 1 — 3 
site (middle) sl z — 6, z' — 4 site (bottom) a z = z = 8 site 



block spins S' are situated on the black circles repre- 
senting the sites of the inflated lattice, while all of the 
nearest neighbors are decimated in the renormalization 
group (RG) transformation. Initially, all spins have the 
same value of spin, so — ^ , so that after one inflation the 

block spin variables are simply = {S^\ , Sjp} = 

{7 sq, 7 Sq, 7 Sq, 7 So, 6so, 4so, 2so}. (Note that for E and F 
sites, the value of z was corrected for the bond dilution). 
In subsequent inflations, one has the following matrix re- 
lation S(") = (4"\ -,4 n) ) = CS^" 1 ), with 
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C = 



(-1 
-1 
-1 
-1 



V o 



8 \ 
8 
8 
8 
2 5 
3 2 
-112 0/ 



(5) 



Note that the number of values of the block spins after 
each inflation does not grow - there are just four pos- 
sible different values of the block spin at any stage of 
inflation. As in the square lattice example, the spins all 
tend to the classical limit as n goes to infinity. In addi- 
tion, the largest eigenvalue of C, 3, is precisely that of 
the square lattice in section IV. 1 ! This eigenvalue, along 
with the corresponding eigenvector gives the flow of ef- 
fective spin values in the limit of large n. Thus for large n 
S (n) ~ 3S(«-!). This is the same spin renormalization as 
that on the square lattice where z is everywhere equal to 
4. In both cases, the spins tend to infinity, i.e. the classi- 
cal limit, under renormalization. On the tiling, moreover, 
the block spins tend to constant relative asymptotic val- 
ues which are site dependent and given by the eigenvector 

fl 1 1 1 1 3 U 



C. Ground state energy of an isolated block 

Consider the configuration of z + 1 spins of Fig^l] in 
which each of the z links represent the same antiferro- 
magnetic coupling J, termed the Hciscnberg Star (HS) 
in |3j . For spin i variables on each site and for a given an- 
tiferromagnetic coupling J between the central spin and 
its z neighbors, the ground state energy can be found 
exactly to be 



,(<>) 



(z) = -J(z + 2)/4 



(6) 



On the octagonal tiling, one has the seven different fam- 
ilies of star clusters on the tiling, with the correspond- 
ing values of z on the right hand side of the equation. 
The superscript "0" indicates that this corresponds to 
the energy of unrenormalized clusters. We also require 
the ground state energy in the case of clusters of spins of 
unequal lengths. The lowest energy state of a cluster in 
which z spins of unequal lengths Si = riiSo are coupled 
with strength J to a central spin Sq — uq-sq is taken to 
be the following generalisation of Eq|S| 



e( J, z, {n}) = -n Q J{^2 n< + 2)/4 



(7) 



i=l 



In the present model, although initially the couplings are 
all equal, after one RG step the couplings take on dif- 
ferent values. Therefore we shall make an approximation 



later that consists of replacing the set of couplings around 
each site by a single locally averaged value. 




FIG. 10: (z+1) spin cluster (Heisenberg Star) 



D. Proliferation of blocks under deflation 

If n = (tia, tib, nc, noi) are the number of blocks in a 
given region of each given type, the number of blocks of 
each type after one deflation is Pn where 



P 



( 1 8 \ 

10 2 5 

10 4 2 

V 1 1 4 / 



(8) 



The largest eigenvalue of the proliferation matrix P 
is equal to 7 so that the total number of blocks 
increases(decreases) with the number m of defla- 
tions (inflations) as 7 m for large m. Notice that the pro- 
liferation of blocks is described by an integer, and not the 
irrational number A 2 w 6.8, each of these numbers being 
the answer to a different question. The former describes 
the rate of growth of a finite system in terms of the num- 
ber of blocks. The latter is the scale factor of the change 
of site density under inflation/deflation for the infinite 
quasicrystal, and this is not restricted to have integer 
values. 



E. Renormalization of links 

There are an infinite number of types of links since each 
link couples two sites that are each unique. However, just 
as we chose to truncate the size of the space of solutions 
by distinguishing only seven types of sites, we can con- 
sider a "minimal" model where it suffices to take into 
account only five kinds of links. These are represented 
in an array j = (j a F 1 j a E 1 jD 1 D 2 , j d 2 f ?j ef) ■ Here, j aF is 
used to denote the link between (A,F), (B,F),(C,F) and 
(£>i,F) pairs. Similarly, j a E denotes the link connecting 
(B,E),(C,E) and (Di,E) pairs. This oversimplification of 
the link classification ignores, in particular, that E and 
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F sites can occur in more than one environment. How- 
ever, in the first approximation, we have assumed here 
that one can treat all the sites of a given family as iden- 
tical out to first neighbors, and this approximation will 
be found post facto to yield reasonably good numerical 
results. 

Note that there arc no bonds linking sites that are 
separated by a distance d s in the original tiling (recall 
that this is the shortest distance possible on the octagonal 
tiling) and the same is true for the sites of the inflated 
tiling since our bond dilution has the effect of decoupling 
such blocks. 

Interblock links are all the links not taken into account 
in the definition of blocks. To find the new effective links, 
one also allows for bond moving, as illustrated by the fol- 
lowing example: consider a central A site surrounded by 
eight D\ -clusters. These transform to an A site with 
eight F sites around it after an inflation. We wish to ob- 
tain the effective link between the central A and one of 
the F sites. The original A site has sixteen links to the 
eight Di-clusters - i.e. it has two links per Di-cluster. 
These two links between the center and each peripheral 
block are of the EF type (see Fig. 8c). Thus the new ef- 
fective coupling between the central A — > A site and each 
of the eight D\ — > F sites around it on the inflated tiling 
is of the aF type. It is antiferromagnetic, like the original 
couplings. One takes into account the spin renormaliza- 
tion factors of the block spins mentioned before, namely 
(;a and £d respectively. The new coupling can then be 
expressed in terms of the previous generation of couplings 
by the equation 



7 -(i) 

JaF 



o,'(0) AO) AO) 

^Jef^a sb 



(9) 



For the second type of links, j aE , one sees that there 
are three EF links joining a4^5 site to a C — + E site, 
so that the new aE link is given by 



7 (1) _ o AO) AO) AO) 
JaE — JEF^A ?C 



The other effective couplings can be written down simi- 
larly, although a problem arises due to the fact already 
mentioned, namely, that E and F sites can occur in more 
than one local environment. Here we chose just one op- 
tion among the several, to write down the new effective 
D 2 F and EF couplings. With this truncation of the link 
relations, we have a system of equations between the five 
old and five new couplings, j*- 1 -* = M^j(°\ where 



M (n) 



/o 






n An) An) 
u SB ?D 

An) An) 








n An) An) 
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°S4 S, 



A sc 
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SB SB 



(ii) 



with the initial condition (taking the zero order coupling 
J=l) = (1,1,1,1,1). 



F. Averaged values of renormalized couplings 

After one inflation, the new tiling has the same ge- 
ometry, with the same relative frequencies of vertices as 
the old tiling, however, the new onsitc spins S' 1 ) and in- 
tersite couplings j' 1 ) are no longer uniform from site to 
site. To proceed, we define averaged quantities - aver- 
aged renormalization factors ^ and averaged couplings, 
for each of the seven types of site. The average couplings 
are easily found, using the local environments listed for 
each of the seven families in Tables 2 and 3. The simplest 
situation occurs for A sites, which have eight A-F links 
surrounding them, so that the average coupling is just 
= j^p . For the six remaining sites we can similarly 
define averaged couplings that are linear combinations of 
the j( n h Dropping the superscripts, we thus have seven 
averaged couplings as follows: 



lA = JaF 

JB = ( 5 iaF + 2j aS )/7 
JC = ( 2 jaF + 3j aE )/5 

Jd 1 = (2joB + .7bb)/3 

Jd 2 = ( 2 jaE + 2.7bb + 2jb 2 f)/5 

j e = {2j aE + 2j EF )/A 

If = U*F + 2j EF )/3 (12) 

Average renormalization factors ^ are analogously 
determined for each of the seven sites, and used to obtain 
the new matrix . This process is repeated and the re- 
sult is a set of recurrence relations = M(")j(") with 
M having the same structure as in Eq.8. One can now 
study the evolution of the matrix M under successive in- 
flations. The maximum eigenvalue of M, 75 w 0.15. This 
results in a power law decay of the couplings for large n, 
since w 75j^™ _1 - ) . The corresponding eigenvector ^5) 



(10) 

determines the fixed point relative couplings. 



G. Hamiltonian of inflated system 

The effective Hamiltonian after a single inflation is now 
written down much as for the case of the square lattice. 
After the first renormalization there arc block spins at 
each of the a-class sites, whose ground state zero or- 
der energies are e- ^ and having new interblock links jt 1 ). 
ffW(iVW,{sf ^{jW}), where ff« has the same form 
as the original Hamiltonian in Eq.2 and = \~ 2 N. 
The original Hamiltonian is thus decomposed into a set 
of independent cluster energies and a set of intercluster 
terms as follows: 



ff = £/;4 0) + " (1) 



(13) 
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where j can take on the values A,B,C or D\, The first 
term is a sum over the energies of Heisenberg stars de- 
fined on the four types of blocks a, given by EqEI or 
equivalently by Eq[7|with = e(l, z, no = 1, Y^ n i = z )- 



V. RESULTS 

We will discuss the calculation of the local order pa- 
rameters and then that of the ground state energy. 



1. Local staggered magnetic moments 

The QMC data in [g give values of local order param- 
eters. These can be defined in terms of the local energies 
around a site i 



-J ^2{Si.S i+s ) 



(14) 



where the sum is over all the nearest neighbors of a given 
site i, and the spin correlations were evaluated in the 
ground state. We have added a factor h per bond (that 
is, the bond energy is shared equally between the two 
sites at each end). The local order parameters are defined 
by 



(15) 



It is the quantity that we now wish to calculate. 

The inflation symmetry of the quasiperiodic system al- 
lows us to define clusters on length scales that increase 
as powers of A 2 . We would like a relation between the lo- 
cal energies ei and the cluster energies, denoted E^(z), 
evaluated as a function of z for bigger and bigger cluster 
size as n increases. The energy per site for a cluster of 
the ith type tends to a certain value in the infinite size 
limit. We propose that this limiting value coincides with 
the local energies calculated by the QMC. This is based 
on the expectation that there is a fixed point distribution 
for cluster energies, like the one found for the block spins, 
and for the averaged couplings. 

The number of terms contributing to the cluster en- 
ergy is governed by the largest eigenvalue of the block 
proliferation matrix P, so that E^ n ' /7™ tends to a limit 
as ti — > oo . It is this quantity that corresponds to the 
numerically evaluated local energies. With this assump- 
tion, the local order parameters at every stage of RG are 
found from 



» _ 



ff(")(z) 
7 n z 



(16) 



We now describe how to calculate the cluster energies 
at each stage of RG. 



Zeroth order calculation 

The zeroth approximation was obtained in Q , the en- 
ergies of the clusters at this order being easily calculated 



using Eq|S| for each of the values of z, 



,(0) 



eC°). The 



values obtained are 



,(°) 

■A > ■ 



(Oh _ rl £ 2 1 1 3 5 
■' 6f >_i 2'4' '4'4'2'4 > [ V 



The staggered moments corresponding to these ener- 
gies are a simple function of z 



i {0) 

loc,s 



(z) 



i(°)(z)/z 



z + 2 
4z 



(18) 



This function is plotted in Figlllb (dashed line). In 
accord with the qualitative trend of the QMC data, it 
shows that toz OCiS decreases with increasing z. With each 
additional bond, the central spin enters into a resonant 
state with more and more neighboring spins, with the re- 
sult that for each individual bond there is less amplitude 
for formation of a singlet. 



A. First order calculation 

The seven averaged couplings at this order have the 
numerical values 



{Ja, -,Jf} = {°- 14 > °- 13 > °- 12 , 0.10, 0.16, 0.24, 0.29} 

(19) 

These averaged couplings are used in the calculation 
of the ground state energy at each of the new clusters. 
This is done using Eq0 along with the block spin values 
for the center and three surrounding blocks deduced from 
Eq.5. The first order Heisenberg star energies for each of 
the seven types of site are thus 



( ® \ 

e (1) 
e (1) 

e Dl 

e (1) 

t D2 

V J 



( e{j A , 8, {no = 7, Ern = 16}) \ 
{n = 7, En, = 18}) 
eGo 6 > {"o = 7, Y,ni = 16}) 
e 0m> 5 ) i n a = 7, Sn, = 14}) 
e (~3D2, 5 > {"o = 6, Em = 19}) 
e(j £ ,4,{n = 4,I]n 4 = 18}) 
V 3, {n = 2, Em = 15}) J 



(20) 



The energy of a cluster at first order, denoted in- 
cludes this Heisenberg star energy and all zero order di- 
agonal terms of the sites belonging to the cluster. These 
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first order energies of the clusters can be expressed as 
follows: 

4> =ef W 0) r , + -yV° 3 ,, (21) 

where j = 1, .., z are the nearest neighbor sites of i, and 
anc(i) denotes the ancestor of site i. This definition takes 
into account the first order star cluster energy for the 
cluster i plus the zero energy term for the center site, 
plus one-half the zero energy terms for the surrounding 
sites. 

To illustrate with an example: consider an A site on 
the inflated tiling, with eight nearest neighbor F sites 
around it. The zero order energy term for an A site is 
the block spin energy of its ancestor A site, namely, ■ 
The zero order energy term for F sites is the energy of 
their ancestor D\ block spins, e^\. Finally, the Heisen- 
berg Star energy for the A site, and with the first order 
effective coupling is . 

Consider another example of an F-site which has three 
neighbors, say an A site and two E sites. The F site 
arises from a D\ site. The zero order block energy asso- 
ciated with it is therefore e^\. Similarly, the ancestors 
of the three neighbors are an A and two C sites. They 
contribute half their block energies, respectively and 
to the total F-cluster energy. The total energy of 
the F-cluster is found by adding four zeroth order terms 
plus the HS energy for F sites, which have a first-order 
coupling 'fp ■ Other cluster energies can be similarly ob- 
tained, and are listed below. 
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obtained from the preceding order zero energy terms by 
use of the proliferation matrix P defined in Eq|S] We give 
the F cluster energy to this order, as an example: 

4 2) = 4 2) + ^ + ^(4 1) + 2 £ g ) ) + 

(^ + 4 0) +84 0) +64°l) (23) 

At third order, proceeding similarly, there will be a 
term in , four terms in e*- 2 -*, and a certain number 
of terms in and e( ). The number of blocks of each 
type can be found easily using the proliferation matrix 
to determine the number of ancestors of each type of 
block. In Fielllh, we have compared the m s obtained af- 
ter zero (the dashed curve) with the results at one and 
two RG steps (open circles and squares). After the sec- 
ond step, the values of m s converge quickly as can be seen 
in Figlllb which shows the third (circles) and fourth or- 
der (squares) results along with the QMC data, 

C. Predictions for the full octagonal tiling 

The limiting values of mi OCjS are clearly below the 
QMC data. This is to be expected, due to the bond 
dilution. One has to correct for the effect of the appre- 
ciable bond dilution occurring at C and D sites in order 
to obtain an estimate of the energy of the undiluted oc- 
tagonal tiling. On the one hand, the bond dilution leads 
to having fewer energy terms in the Hamiltonian and con- 
sequently underestimating the cluster energies. On the 
other hand, the loss of bonds is partly offset by the fact 
that the dilution tends also to suppress frustration and 
raise the local order parameter. An ad- hoc way to put 
back the "missing bond-energies" is to add in half of the 
missing link energies at each of the C and D sites. This 
is easily done here by adjusting the z values at each of 
the sites, zc goes up from 5 to 5.5 while zdi is increased 
from 3 to 4. Using this ad-hoc procedure we can get 
estimates for m s values on the original octagonal tiling. 
The grey squares of Fig lllb were obtained by adjusting 
the n = 4 data in this way. As the figure shows, this 
procedure yields a fairly good agreement with the QMC 
data. The same procedure is used to obtain the ground 
state energy estimate of the full octagonal tiling in the 
next section. 



B. Second order calculation and higher orders 

For n — 2, the energies of the seven clusters for the 
twice-inflated tiling can be written out in terms of the 
energies e^ k \z) (k — 0,1,2). It is easy to obtain the 
explicit expressions since it suffices to increase all the 
superscripts in Eq|22by one (so for example the e^' be- 
come e\ )■ The zero order energy terms are also easily 



2. Ground state energy 

The ground state energy Eg is the sum over all blocks 
at all orders, of the block energies. At zero order the 
number of blocks of z-spins iV^ ' (z) = Nfi ( i.e. propor- 
tional to the original frequencies of occurrence given in 
Eq.l). The density of vertices decreases with each infla- 
tion as A 2 , so that 
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FIG. 11: m s (z) values versus z obtained for increasing or- 
ders of RG. The zero order analytical curve is indicated by 
a dashed line in each figure, (a) 1st (circles)and 2nd (rect- 
angles) order RG. (b) 3rd (circles) and 4th (rectangles) order 
RG and QMC data for the full octagonal tiling, (c) Adjusted 
4th order data (grey rectangles) and QMC data (grey circles). 



E /N = J2fi(< 



x 2 



' A 2 ™ e% 



(n) 



.) (24) 



The block energies e^> are the energies of blocks with 
a spin Sq 71 ^ at the center, with effective couplings j*-™) 



to the S. 
gives e 



(n) 



surrounding spins. The series for the energy 
—0.51. We can estimate the effect of bond 



dilution, as was done for the local order parameters. Us- 
ing the corrected values of z explained in the last section, 
one finds an adjusted ground state energy of about —0.59. 
This value of the GS energy is significantly smaller in ab- 
solute value than the value deduced from the QMC data 
in . We recall that this was true of the square lattice 
calculation as well. In that case, the RG calculation of 
Delgado and Sierra was already noted in [lfj to underes- 
timate the bonding energies of pairs of spins because of 
the inadequacy of first order perturbation theory around 
the Neel state. The same is presumably true of our RG 
on the octagonal tiling. For the former case the RG cal- 
culation was compared with the terms of a 1/5* expansion 
of the ground state energy, and shown to lack the sub- 
leading order term, resulting in the observed discrepancy 
of values. 

On the square lattice, eo has been determined numeri- 
cally j2{| to high precision to be —0.6694, while finite size 
scaling for the tiling [29J obtains a value of —0.6581. The 
closeness of the values obtained for these two very differ- 
ent problems is rather surprising. It is a probable that 
this close proximity of values is due to the fact that the 
octagonal tiling, with its two sublattice structure and its 
average coordination number of 4. The differences must 
arise from the next nearest neighbor distributions which 
differ for the two systems, although this remains to be 
verified by explicit calculation. 



VI. DISCUSSION AND CONCLUSIONS 

In conclusion, we have presented an approximate RG 
scheme for ground state properties of a two-dimensional 
quasiperiodic tiling that can be solved after bond dilu- 
tion. Other approximations involve the truncation of the 
number of distinct sites and the number of distinct links, 
and replacing local couplings around sites by average val- 
ues in order to simplify the effective Hamiltonian after ev- 
ery inflation. The results obtained for the diluted tiling 
were used to get estimates for the undiluted tiling. De- 
spite these approximations, we believe the model solved is 
close to the perfect two dimensional quasiperiodic struc- 
ture, and it allows for a rather detailed solution of real 
space properties of these heirarchical structures. The 
results obtained by RG for local order parameters are 
close to those calculated for the full undiluted model, af- 
ter our adjustment procedure. It thus appears that the 
model takes into account the most relevant aspects of the 
quasiperiodic geometry of the octagonal tiling. 

The RG method presented is less good at obtain- 
ing the ground state energy, similar to the situation al- 
ready noted for the square lattice by Sierra and Martin- 
Delgado, who showed that a better result is obtained by 
going to second order of perturbation theory to obtain 
the effective Hamiltonian after renormalization. Con- 
cerning the proximity of values of the ground state en- 
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ergy in these two systems, our calculation is not accu- 
rate enough to explain this observation. A calculation to 
higher order would involve further nearest neighbor sites, 
improve the energy estimate and perhaps help to explain 
the small energy difference between the tiling and the 
square lattice. It would be interesting as well to com- 
pare results for other bipartite two dimensional tilings, 
including the Penrose tiling. 

The zero temperature magnetic state of this quasiperi- 
odic Heisenberg antiferromagnet has a structure factor 
with peaks that can be indexed using the four dimen- 
sional indexing scheme (see Appendix). The positions 
of the peaks is very simply related to the positions of 
the peaks of the paramagnetic state: they are situated 
halfway in between. In other words, the paramagnet is 
indexed by four integers, while the antiferromagnet has 
half-integer entries, corresponding to the antiferromag- 
netic vector q = {^,5,^,5}- This is the quasiperiodic 
analogue of the square lattice where just such a shift 
occurs in reciprocal space and corresponds to the anti- 
ferromagnetic vector q = {^,5} (see [3tj for a discus- 
sion along with a simple one dimensional version of a 
quasiperiodic antiferromagnet). The real life quasiperi- 
odic compound ZnMgHo was studied by neutron scat- 
tering and shown to have short range antifcrromagnetic 
correlations below about 20K. These correlations lead to 
a magnetic superstructure that is, as for our two dimen- 
sional model, shifted with respect to the paramagnetic 
state. The antiferromagnetic vector that best fits the 
data has a more complicated value than the simplest form 
for a 3d quasiperiodic antiferromagnet (c^ = i, i = 1, 6). 
This is because the magnetic unit cell is much larger for 
the three-composant system, due to the fact that only the 
Ho sites carry a magnetic moment, resulting in smaller 
spacings between peaks in reciprocal space. 

Finally, the RG scheme presented here can be adapted 
to discuss other discrete quasiperiodic models, such as 
tight-binding models for electrons hopping between ver- 
tices of the tiling. It should provide a useful theoretical 
framework for describing quasiperiodic tilings in general. 

APPENDIX. THE CUT-AND-PROJECT 
METHOD. 

1. One dimensional example 

The cut-and-project method of obtaining quasiperiodic 
tilings is easiest to illustrate in the case of the celebrated 
one-dimensional tiling - the Fibonacci chain(see Luck's 
review in [22(). The Fibonacci chain comprises two basic 
tiles or line segments of two different lengths, "long"(L) 
and "short" (S) arranged in a deterministic sequence. The 
Fibonacci chain can be generated iteratively from a sin- 
gle S segment using the following substitution rules: re- 
place each S by an L, and each L by SL. Two succes- 



sive segments of the infinite chain are shown below to 
illustrate the substitution rules (dashed lines represent 
L, thick lines S) 



FIG. 12: two segments of the Fibonacci chain 



The Fibonacci sequence of segments, or tiles, can 
be generated by projecting selected edges of a two- 
dimensional square lattice onto the one dimensional 
"physical space" E± as shown in Fig El The vertical 
and horizontal edges project onto the S and the L tiles 
respectively. The orientation of E\ is given by tan -1 1/r 
(where r = (V5+ l)/2 is the golden mean, a solution 
of r 2 — t — 1 = 0), an irrational slope, so the tile se- 
quence never repeats. The edges selected for projection 
onto E\ obey the following condition: the projection of 
the edge onto the perpendicular space E2 must fall within 
the "window of selection" W (indicated by the thick line 
segment representing the projection of the unit square 
shown in grey) . A finite sequence of twelve edges satisfy- 
ing this condition are shown in bold in Fig ll3l and they 
result in the projected structure ...LSLLSLLSLSLL... 




FIG. 13: Cut-and-project method. Selected edges of a square 
lattice are projected onto the parallel space (-Ei). The edges 
that are selected have perpendicular space (-E2) projections 
that fall within the segment marked W. 



2. Two dimensional case 

In analogy with the one dimensional case, the octag- 
onal tiling is obtained from the projection onto E\\ (the 
physical two-dimensional space) of a subset of vertices of 
a four-dimensional cubic lattice. The subspaces Ex,E% 
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are now two-dimensional, and are invariant under eight- 
fold rotations in the four dimensional space. The orien- 
tation of the physical plane Eh is given by the number 
A = 1 + y/2, one of the solutions of A 2 - 2A - 1 = 0. The 
tiles, which are projections in this plane of the 8 faces of 
the 4d cube, are squares and 45° rhombuses. The ver- 
tices (and edges) that are selected for projection in the 
two dimensional perpendicular space E2 must fall within 
the window of selection shown in FigEl This octagon- 
shaped area is delimited by the projection of the sides 
of a four dimensional unit cube. The octagonal tiling 
has by construction the eight-fold symmetry in the weak 
sense already described. There are six different kinds of 
nearest neighbor configurations for vertices of the tiling, 
denoted A through F as shown in Fig. 2. 

A. Domains of acceptance 

The maximum coordination number z — 8 corresponds 
to the eight-fold symmetric A sites. These sites possess 
perpendicular space projections that always fall within 
the central octagon labeled A in Fig^] Similarly, the 
projection of all z = 7 sites falls within one of the eight 
triangular regions adjoining the central octagonal do- 
main. The remaining z values likewise correspond to the 
domains labeled accordingly in Fig^| which has exact 
eight-fold symmetry The ratio of side lengths of similar 
polygons are either A, or A 2 . 









F \ 










E / 


/ c 


/ \ 

A 

















FIG. 14: Projection into perpendicular space of vertices of 
the octagonal tiling. Domains corresponding to the six fam- 
ilies are labeled ( eight-fold symmetry determines labels of 
unmarked domains) 



B. Inflations and deflations 

For the octagonal tiling, the inflation transformation is 
given by a 4 x 4 matrix acting in the four dimensional cu- 



bic lattice, satisfying U 2 — 2U —1 = and having entries 
of or 1 only. Projecting the subset of points selected by 
11 leads to a bigger tiling of the same type as the original 
one. Only the highest z sites remain selected, while the 
others disappear. The sites that remain are those within 
the middle octagon, the a family: A, B, C and D\. The 
sites that disappear correspond to the region outside the 
middle octagon. The perpendicular space representation 
also allows to determine rapidly the new z values of the 
sites that remain: one simply redraws the acceptance do- 
mains of Fig^Jafter rescaling, inside the middle octagon. 
Thus a point that was previously in the D\ domain will 
find itself in the F domain, and the C domains map into 
E domains. A sites remain A sites if they are close to the 
center of the diagram, otherwise they become one of the 
other a sites after inflation. The four categories of A sites 
mentioned in section II. 3 differ by their distance from the 
origin in perpendicular space. The perpendicular space 
projection of a site determines its evolution under infla- 
tion - the closer a site is to the center of the octagonal 
selection window, the longer it remains an A site under 
successive inflations. Also, one clearly sees that D-sites 
come in two types, with different perpendicular space do- 
mains. The following table resumes the old and new site 
types after inflation: 

A — > A or B or C or Di 

B — > Da 

C -> E 

Di -> F 

The number of sites per unit area is reduced by the 
scale factor A 2 after each inflation, and the relative fre- 
quencies of occurrences of each of the seven families of 
sites is invariant. The frequency of occurrence of the ith 
family is proportional to the area occupied by that fam- 
ily in the perpendicular space projection. It can thus 
be easily verified using Fig^^that these frequencies are: 
Ia = A" 4 ;/* = A- 5 ;/ c = 2A- 4 ;/ m = A~ 3 =; f E = 
2A~ 2 ;/p = A~ 2 . The average coordination number is 
exactly 4, as can be checked using the frequencies given. 
The interested reader can find these and other impor- 
tant geometric and algebraic properties of the system 
described for example in |23|. 

C. Reciprocal space and structure factor 

The diffraction peaks of the octagonal tiling are found 
at positions given by projections into 2d of reciprocal lat- 
tice vectors a.; of the 4d cubic structure. The intensities 
of the peaks are not uniform however, but depend on the 
Fourier Transform (FT) of the finite selection window. 
The main features of the diffraction pattern are thus (see 
Belin et al in 2] for more on the topic) 

• The peaks have an eight-fold symmetry around the 
peak at the origin. 
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• peaks occur at positions corresponding to the set of 
integers h, k, to, n representing the projection into 
the 2d plane of the 4d vector q = h&i + fca 2 +ma 3 + 
rta 4 . That is, one can index peaks by a set of four 
integers. 

• Intensities are highly dependent on the value of q 
since the FT of the selction window is oscillatory 
and long ranged. The set of eight most intense 
peaks nearest the origin is used to define a quasi 
Brillouin zone for the tiling. 

D. Approximants and some of their properties 

Numerical studies of quasiperiodic systems are per- 
formed on finite pieces of the infinite system. In par- 
ticular, it has been pointed out that periodic boundary 
conditions are preferable to open or closed boundary con- 
ditions in terms of eliminating spurious states and eigen- 
values. A periodic approximants is a structure that can 
be periodically continued and can be augmented in size 
so as approach arbitrarily close to the perfect infinite 
structure. 

This is, again, easiest illustrated by going back to the 
Fibonacci chain. In the cut and project technique, it 
should be clear that if one tilts the irrationally oriented 
selection strip away from the special angle, one will ob- 
tain a periodically repeating chain every time the slope 
is rational. 

r _1 has a scries of approximants given in terms 
of the Fibonacci numbers as follows: {a i, ...} = 
{1, |, |, |, -p^-...}, where Fk is the fcth term in the 
Fibonacci sequence defined by the recurrence relation 

Fk+i = Fk + Ff.-i 

F =F 1 = 1 (25) 

with Fq = F i = 1. By increasing the value of the 
denominator of the rational number F n /F n+ i - i.e. by 
choosing increasingly longer approximants of the golden 
mean - one will get a structure of period F n+ 2- The 
finite sequences of L and S within the approximant are 
the same as those found in an infinitely long chain. 

For the two dimensional case, Ref.j2jJ describes how 
to obtain square approximants to the octagonal tiling 
by the projection method. These are obtained from the 
approximants to the silver mean which depend on ratios 
of the so-called Octonacci sequence : Ak = Ok+i/Ot 
where 

O fc+ i = 20 k + O fe _i 
Ol = 1; 2 = 2 (26) 

with Oi = 1;C>2 = 2. These are the finite size sys- 
tems used for a number of numerical studies including 



the quantum Monte Carlo calculations. The first few 
square approximants have the following sizes 



k 


2 


3 


4 


5 


N k 


239 


1393 


8119 


47321 



Table A.l. Number of sites in the first four square 
approximants 

We will list some features of these approximants that 
may be important to bear in mind depending on the mod- 
els studied. 

• Reflection symmetry (exact) with respect to the 
bottom left-top right diagonal 

• 90° Rotation symmetry around the center (approx- 
imate) 

• Odd parity of repetition. By this is meant that one 
changes sublattice when one goes from a site to its 
first periodic repetition along either x or y direc- 
tions. For a number of numerical calculations it is 
easiest to restore the bipartite property by taking 
a system size doubled along both directions (i.e. 
quadrupled unit cell with respect to the sizes given 
in the Table). 

• Inflation relation between approximants. Fig 1 151 
shows a small approximant superimposed on the 
next largest one. 




FIG. 15: Superposition of two successive approximants 
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